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Abstract 

PHENIX has measured the centrality dependence of charged hadron pT spectra from Au+Au 
collisions at \/{snn) = 130 GeV. The truncated mean px decreases with centraUty for 
Pt > 2 GeV/c, indicating an apparent reduction of the contribution from hard scattering 
to high pt hadron production. For central collisions the yield at high px is shown to be 
suppressed compared to binary nucleon-nucleon collision scaling of p+p data. This sup- 
pression is monotonically increasing with centrality, but most of the change occurs below 
30% centrality, i.e. for collisions with less than ~140 participating nucleons. The observed 
Pt and centrality dependence is consistent with the particle production predicted by models 
including hard scattering and subsequent energy loss of the scattered partons in the dense 
matter created in the collisions. 



1 Introduction 



Particle production at large transverse momentum (p^) provides a new tool to study 
hot and dense nuclear matter created in high energy nuclear collisions. In nucleon- 
nucleon collisions, hadrons with pt > 2 GeV/c are believed to originate mostly 
from the jet fragmentation of constituent partons, quarks and gluons, that were 
scattered with large momentum transfer [1]. In nuclear collision these hard scat- 
tering processes between constituent partons occur early compared to the lifetime 
of the strongly interacting matter. Thus the hard scattered partons may traverse the 
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highest energy density matter produced. Theoretical studies of the propagation of 
partons in high density matter suggest that they lose a significant fraction of their 
energy through gluon bremsstrahlung [2] and that the energy lost reflects the den- 
sity of color charges in the matter through which they pass [3]. The energy loss 
reduces the momenta of the partons, which results in a corresponding reduction of 
the momenta of the fragmentation products [4], observable as a reduced yield of 
high pt hadrons. 

The first measurements of hadron spectra at the Brookhaven National Laboratory 
Relativistic Heavy Ion Collider (RHIC) facility indicate a suppression of high-pr 
hadron production in central Au+Au collisions relative to a binary collision scaling 
of p+p and p +p data [5,6]. No suppression is found for peripheral Au+Au col- 
lisions. So far no unique explanation of this apparent absence of the expected jet 
contribution to the pt spectrum above 2 GeV/c has been identified. 

Models of parton energy loss can reproduce the observed suppression in central 
Au-i-Au collisions [7,8,9]. Other final state effects such as re-scattering of hadrons 
originally produced via the jet fragmentation have been proposed to explain the 
suppression [10]. It should be noted that models invoking thermal hadron produc- 
tion combined with collective transverse expansion of the reaction volume success- 
fully describe the transverse momentum distributions of identified hadrons up to 3 
GeV/c [11,12]. However, the mechanism of equilibration, which requires a reduc- 
tion of the high pt particle yield, is not specified in these models. 

Alternatively, the initial state may be modified such that the number of hard scat- 
terings is reduced. It is well known that nuclear modifications of the parton dis- 
tributions exist [13]. These modifications cannot explain the suppression, since in 
the kinematic range of the measurements anti-shadowing enhances the parton dis- 
tributions in nuclei [14,15,16]. However, models using a classical QCD picture of 
a highly relativistic nucleus [17,18] suggest that gluon distributions are saturated 
for momenta below a scale Qs and thus reduced compared to expectations based 
on perturbative QCD [19]. As a consequence, a considerable suppression of hadron 
production might be expected even well above Qs [20]. 

In this letter we present the centrality dependence of the suppression of the high-pT 
hadron yield to provide new experimental constraints on theoretical descriptions. 
These data are complementary to the previous study of the absolute yields [5] and 
have different systematic errors. 



2 Experimental setup and data analysis 

The results are obtained from 1.4x 10^ minimum bias Au+Au collisions at y/s^=130 
GeV recorded by the PHENIX experiment during the Run-1 operation of RHIC 
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(August - September 2000). Details on the PHENIX detector and its configuration 
in Run-1 operation can be found in [5,21]. 

In PHENIX semi-inclusive charged hadron spectra are measured over the range 
0.5 < pt < 5.0 GeV/c in the east central arm spectrometer using data from a drift 
chamber (DC) and two segmented cathode pad chambers (PCI and PC3), located 
outside of an axial magnetic field at a radial distance of 2.2, 2.5 and 5 m from 
the beam axis. The detectors cover an azimuthal acceptance of 90° and a pseudo- 
rapidity range of \r]\ < 0.35. In this analysis an additional fiducial cut I77I < 0.18 is 
applied to guarantee homogeneous track acceptance for collisions within \zi,tx\ < 
30 cm of the nominal interaction point. About ~25% of the azimuthal acceptance 
is covered by a time-of-flight system which allows proton identification out to 3.5 
GeV/c, where the measurement is limited by statistics [22]. 

A pair of beam-beam counters provides the vertex position along the beam di- 
rection (z). Each charged track is reconstructed from the DC measurements of its 
projection into the bend plane of the magnetic field and two space points provided 
by PCI and PCS. The unphysical background, resulting from false associations of 
drift chamber projections with pad chamber points, is estimated and subtracted by 
forming artificial events with the locations of pad chamber points inverted around 
the symmetry axis of the spectrometer. Physical background from decays in flight 
and photon conversions close to the DC, which only partially traverse the field and 
thus mimic high momentum tracks, are removed by requiring the track to point 
back to the event vertex within \zytx\ < 2.5 cm. The remaining background level 
is negligible below 4 GeV/c and less than 40% at 5 GeV/c; this upper estimate is 
included in the systematic errors. 

Corrections of the data are determined by tracing individual particles through a full 
GEANT simulation, simulating the detector response, then merging this response 
with that of all particles from a real event and passing the composite event through 
the PHENIX reconstruction software. The average track reconstruction efficiency 
in the active detector area is larger than 98% in peripheral collisions and decreases 
to 68±6% for central collisions. The corresponding correction is shown on the up- 
per left hand side of Fig. 1. The full correction also depends on px- It is plotted 
for peripheral collision on the upper right hand side of the figure. Between 0.8 and 
2.5 GeV/c the correction factor varies slowly with p^. Its value of ~ 25 corrects for 
geometrical acceptance (A0 = 7r/4; Ar] = 0.36), dead areas of the detectors ( 45% 
DC, 5% PC's), and losses due to 2a track matching cuts. At lower pt the correc- 
tion increases reflecting the gradual loss of acceptance. At higher pr the observed 
particle yield is artificially increased because of the finite momentum resolution 
(Sp/p ~ 0.6% ® 3.6% p (GeV/c)) and therefore the correction function decreases. 
Since this correction depends on the spectral shape of the true pt distribution it was 
determined iteratively. At 5 GeV/c the correction is reduced by a factor of ~2. The 
systematic uncertainties are indicate by the dashed lines; they are also tabulated in 
Table 1. As shown in the lower part of Fig. 1 correction factorizes into functions of 
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centrality (i.e. detector occupancy) andpr within 2% systematic uncertainty in the 
range from 2 to 5 GeV/c. 

Events are selected according to centrality following the procedure described in [5]. 
Six exclusive centrality bins are established using the energy measured in two zero- 
degree calorimeters and the number of charged-particles detected in the two beam- 
beam counters. A Monte-Carlo simulation using measured nucleon density distri- 
butions calculated in the Glauber eikonal approximation was used to estimate the 
average number of binary nucleon-nucleon collisions (Ncou) and the correspond- 
ing average number of participating nucleons (Npart) for each bin. The results are 
quoted in Table 2. 



3 Results 



Fig. 2 presents charged hadronpr spectra for the six centrality bins. For peripheral 
collisions the spectra are more concave than those for central collisions. This shape 
difference is seen more clearly by taking the ratio of the spectrum for each centrality 
bin to the minimum-bias spectrum, as shown on the right hand side of Fig. 2. In 
these ratios most systematic errors cancel or affect the overall scale only. The ratios 
for the central bins are almost independent of px since central collisions dominate 
the minimum bias yields. The peripheral bins show a decreasing ratio between 0.5 
and 1.5 GeV/c thus in peripheral collisions the yield in this region decreases more 
rapidly with increasing px than in central collisions. For pt above 1.5 GeV/c this 
trend is inverted. 

Before analyzing the centrality dependence in more detail we demonstrate that the 
observed suppression of the yield at high px does not result from a reduced yield of 
protons and anti-protons [22]. To evaluate the effect of the (anti-)protons, we plot 
in the top panel of Fig. 3 the pr dependence of p/h, the ratio of proton plus anti- 
proton yields to the total charged hadron yield for minimum bias collisions, which 
increases steadily. Above 1.5 GeV/c the ratio seem to saturate reaching a value of 
~ 0.5 arround 3 GeV/c. In the bottom panel of Fig. 3 we show p/h for p^ above 1.8 
GeV/c as a function of centrality. Since there is clearly no significant decrease of the 
p/h ratio with either centrality or pr, the observed hadron suppression is not due to 
a larger suppression of the (anti-)proton component than that of the mesons. On the 
contrary, within the range of the presented measurement the apparent slight increase 
m p/h with Npart could indicate a larger suppression of the meson component 
relative to all charged hadrons. 

To evaluate the change of the hadron spectra more quantitatively we calculate the 
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truncated average pr'- 

fjj^in Pt ■ dNl dpT 

/ trunc\ — fx ' ^mm /i\ 



for p™*" = 0.5 GeV/c and for p™" = 2.0 GeV/c for each central! ty selection ^ . In 
Fig. 4 (p^""'^) is plotted as a function of A'pa^t- The value of (p^""'^) is insensitive 
to the normalization of the spectra. Systematic uncertainties of ~20(50) MeV/c for 
p™*"^ = 0.5(2.0) GeV/c result from an 1-2% uncertainty on the momentum scale 
and the Pt dependent uncertainties on the correction function. Since the centrality 
and Pt dependence of the correction factorize the error on (p^""*^) is independent 
of centrality to better than 2 MeV/c. 

The for pV^"^ = 0.5 GeV/c increases with Npart similar to the average pt 

of identified charged hadrons [22]. For higher p™" = 2.0 GeV/c the (p^^""'^) drops 
by ~60 MeV/c with Npart, distinctly different from the expected increased role 
of hard particle production in the more central collisions. Absent any collective 
effects, the hard scattering contribution should increase relative to soft production 
by the factor Ncoii/ Npart which grows from 1.5 to ~ 6 from peripheral to central 
collisions. Since the relative contribution of the hard component to the spectrum 
increases with px, this should lead to a rise of (p^""'') for sufficiently large p^^"". 
The drop of {pr) therefore indicates the suppression of the high pr relative to the 
low Pt hadron yield independent of systematic errors associated with the absolute 
normalization of the spectra or any nucleon-nucleon reference distribution. 

Changes of the hadron spectra at high pt are often presented in terms of the nuclear 
modification factor Raa- This measure relies on the absolute normalization and 
a reference and therefore has intrinsically larger systematic uncertainties, but it 
allows to quantify the suppression. We have calculated Raa for each centrality bin 
as: 



For the N+N charged hadron cross section we use a power-law parameterization 
l/vr d'^a/dpl = A/(l + Pt/Po)", with A =330 mb/(GeV/c)2, p^ = 1.72GcV/c, 
and n = 12 A. The parameters were obtained by interpolating p+p and p+p data to 
y/s — 130 GeV as described in [5]. In Fig. 5 the Raa{pt) values for all centrality 
bins excluding the most peripheral one are compared to the central (0-5%) bin. The 
systematic uncertainties in the normalization of the data, in Ncou, and in the N+N 

^ The value of (pif^^'^) is closely related to the local inverse slope which is slightly smaller. 
The conversion to local slope depends on the spectral shape and also on p™"*. For an 
exponential spectrum with an inverse slope of 350 MeV/c the conversion is approximately 
-80(-60) MeV/c forp^^" = 0.5(2.0) GeV/c. 



8 



reference (20%) result in overall systematic errors of about 41%, 34%, 31%, 31% 
and 30% for centrality bins 2-6, respectively. The errors are quoted for the range 
from 1 to 3.5 GeV/c; they increase somewhat towards higher p^. Comparing Raa 
for a given centrality bin to the most central bin, the systematic errors reduce to 
~ 6% for bin 5 (5-15%) and nearly 27% for bin 2(60-80%). They are dominated 
the uncertainty in Ncoii listed in Table 2. 

For the 60-80% centrality bin, Raa increases with px and reaches unity at high p^. 
In comparison to the 60-80% bin, the Raa values for the most central bin remain 
significantly below unity at a value of 0.55 for px > 2 GeV/c. At high pt approx- 
imately constant Raa values are detected in all centrality bins. The high pr values 
decrease monotonically with centrality, clearly indicating that the magnitude of the 
suppression of high px hadrons increases with centrality. This is shown clearly in 
the upper part of Fig. 6, which presents Raa obtained for the three pt bins 1.6 to 
2.6 GeV/c, 2.6 to 3.6 GeV/c, and above 3.6 GeV/c as a function of centrality. For 
central collisions we observe a suppression of about a factor of 2 ± 0.6 compared 
to binary collision scaling. Relative to peripheral collisions the suppression factor 
in central collisions increases withp^ from 1.25 ± 0.2 to 1.5 ± 0.2 to 1.8 ± 0.3 for 
the three pt bins, respectively. 

We note that for peripheral collisions the data do not indicate a significant increase 
of Raa above unity, unlike data at lower energies [23]. However, such an increase, 
attributed to initial state scattering, the Cronin effect [24], may well be consistent 
with the peripheral data due to the large systematic uncertainty of the Raa scale. 
While the relative difference between the peripheral and central spectra increases 
with Pt, the roughly constant nuclear modification factor at large pr suggests an 
approximately pr independent suppression of hard scattering contributions over the 
range 2 < pr < 4.5 GeV/c. 

The physics that controls the production of high-p7- particles or the suppression of 
the hard scattering yields in the measured px range may not depend directly on 
Ncoii- Thus, we have calculated a different ratio, Ra'a' defined similarly to Raa 
but with Ncou replaced by the number of participant pairs, Npart/2. If particle pro- 
duction increases proportional to the number of participants, R^aa — 1- 

The obtained R^aa values are shown in Fig. 6 (bottom) for the three pt bins used 
above. The values of R^aa larger than Raa by a factor equal to 2Ncou/Npart 
the average number of nucleon-nucleon collisions suffered by each participant. For 
all Pt bins the yield per participant is consistent with unity for peripheral collisions 
as expected since peripheral collisions should closely resemble N+N collisions. 
For central collisions R^aa increases to approximately three. Most of this change 
occurs in the range of Npart from 40 to 140. For larger Npart the yield in the highest 
Pt bin is approximately constant while in both lower pt bins it increases by 20 to 
30%. 
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4 Concluding discussion 



In this paper we have presented the centraUty dependence of charged hadron pt 
spectra focusing on the behavior of the spectra at high pt. A striking change of 
the spectral shape is observed when comparing spectra from different centrality 
selections. For peripheral collisions the spectrum exhibits a pronounced concave 
shape which is modified towards a more exponential spectrum as the centrality 
increases. The observed lack of variation with centrality in the proton to charged 
ratio at large transverse momenta indicates that the modification is not due to a 
change in the relative yields of protons. 

We observe a decrease of (p^'"^'^) for px > 2 GeV/c with increasing centrality, 
which is distinctly different from the increase of (pr) and demonstrates the suppres- 
sion of the high px hadron yield independent of systematic errors associated with 
the absolute normalization of the spectra. The data are not consistent with binary 
collision scaling of hard scattering processes, which would results in an increase 
of (p^''"'^) . If the Pt spectra above 2 GeV/c are strongly affected by collective mo- 
tion of matter before freeze-out, we would also expect an increase (p^""'^) since 
the corresponding flow velocities should increase in more central collisions [25]. 
Similarly, if gluon saturation is important for particle production in the pr range 
above 2 GeV/c, {ptj:'^"-'^) should increase with increasing Npart due to the predicted 
logarithmic increase of Qs [26]. In contrast, the data are consistent with models 
assuming energy loss of hard scattered partons, which results in an increasing re- 
duction of the hard scattering contribution to the hadron spectrum with increasing 
centrality of the collisions [7,8,9]. It remains to be seen whether this explanation is 
unique. 

Comparing the measured differential yields in five centrality bins to an Ncou scaling 
of the N-i-N reference yields we see a suppression of the yields in central collisions 
at high Pt, consistent with the results in [5,6]. In the 0-30% centrality range (bins 
1-3) the suppression is approximately independent of pr for 2 < pt < 5 GeV/c 
at a value of Raa ~ 0.6 and simultaneously nearly independent within 20% of 
centrality. The suppression sets in gradually with the largest change occurring over 
the 30-60% centrality range. This centrality bin covers a broad range of collision 
geometries. Whether the change is continuous or exhibits a threshold behavior, as 
predicted in [27], cannot be judged from the present data. The observed suppression 
is consistent with parton energy loss scenarios. In these models, the value of Raa 
and its pt dependence are very sensitive to the actual energy loss prescription. Due 
to the large systematic errors on the Raa scale, the contribution from the protons 
and the limited pt reach of the data presented here, we can not distinguish between 
the different energy loss prescriptions on the basis of Raa- 

In summary, a detailed analysis of the centrality dependence of charged particle 
data from Au-Au collisions at ^/snn — 130 GeV measured by PHENIX reveals 
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interesting features of the observed high pt hadron suppression. The decrease of the 
average px with increasing centrality seems to favor models of particle production 
that consider energy loss effects, rather than saturation- or hydrodynamics-based 
approaches for this px range. The suppression sets in gradually with the largest 
changes occurring for peripheral collisions with less than about 140 participating 
nucleons. From there on it does not change substantially towards more central col- 
lisions. 
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Note Added in Proof 

After submission of our manuscript, data from Au-Au collisions at ^/sNN = 200 GeV 
were presented by the PHOBOS collaboration[28,29], with an emphasis on the lack 
of variation in the scaled yields with Npart for Npart > 65. Keeping in mind that 
results from different energies are not directly comparable, we note that the broader 
range in Npart presented here shows that this effect does not apply over the entire 
range of centralities, and that when normalized with the appropriate p-p yields, is 
simply an aspect of the smooth variation of -R^a* visible in Fig. 6. 
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. correction (peripheral)/correction (central) 



0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 5.5 

Pj (GeV/c) 



Fig. 1. Functions used to correct the charged particle pT spectra. Upper left panel shows 
the centrality dependent correction c{Npart) and the right panel shows the pT dependent 
correction c{pt)- The systematic uncertainties are indicated by the dashed lines. The two 
corrections factorize, so that for any centrality the full correction function is given by 
c{pt) X c{Npart)- The accuracy of this factorization is demonstrated in the lower panel. 
The ratio of the full correction for central collisions (top 5%) to the correction for single 
particle events varies by less than 2% above 1 GeV/c. 
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Fig. 2. The left panel shows pT spectra of charged hadrons from six Au+Au centrality 
selections. Error bars indicate statistical errors only. The pT dependent systematic errors are 
independent of centrality and not shown, they are given in Table 1. The centrality dependent 
errors are less than 10% and small compared to the symbol size. The right panel shows the 
ratio of each of the centrality selected px spectra to the minimum bias spectrum. Ratios for 
peripheral selections are scaled for clarity. Dotted lines indicate the average ratios for each 
centrality selection. 
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Fig. 3. The ratio p/h represents tlie proton plus anti-proton yield relative to the total 
charged hadron multiplicity. The top panel shows the px dependence of p/h for minimum 
bias events. In the bottom panel we show the centrality dependence of p/h for pT > 1-8 
GeV/c. Only statistical errors are shown. 
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Fig. 4. Centrality dependence of (p^?""'^), the average of charged particles with px 
above a threshold p™*" minus the threshold p™*". Shown ai^e values for two p™*" cuts, 
one at px > 0.5 GeV/c representing all data presented in Fig. 2 and the other one at 
Pt > 2 GeV/c. Only statistical errors are shown; see the text following Eq. 1 for a dis- 
cussion of the sytematic errors. 
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Fig. 5. Nuclear modification factor (Raa) for the 60-80%, 30-60%, 15-30%, and 5-15% 
centrality selections compared to the one for the most central sample (0-5%). Due to in- 
sufficient statistics Raa is not shown for the 80-92% sample. The solid error bars on each 
data point are statistical. The systematic error between the more peripheral and the central 
sample are given as brackets for the more peripheral data points. The error bar on the left 
hand side of each panel indicates the overall systematic error on the Raa scale. 
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Fig. 6. The top panel gives the nuclear modification factor Raa for three exclusive pT 
regions as a function of the centrality of the collision. The lower panel shows essentially the 
same quantity but normalized to the number of participant pairs rather than to the number 
of binary collisions. The dotted line indicates the expectation for scaling with the number 
of binary collisions (top) or with the number of participants (bottom). Only statistical errors 
are shown. The systematic error on the scale and the centrality dependence are identical to 
the errors shown in Fig. 5. These errors are correlated, i.e. take their maximum or minimum 
value simultaneously for all centrality and pT selections. In addition, there are also px 
dependent systematic errors, which are given in Table 1 . The systematic errors do not alter 
the trends in the data. 
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PT (GeV/c) Strack (%) ^decay (%) ^reso (%) Sbgr (%) total 



1 


±13.5 


+ 10 


±0 





-13.5 +16.4 


2 


±13.5 


+ 5 


± 1 





-13.7 +14.4 


3 


±13.5 


+ 2.5 


±4 


-1.6 


-14.2 +14.2 


4 


±13.5 


+ 1.25 


±9 


-11.5 


-20 +16 


5 


±13.5 


+ 0.6 


± 15 


-40 


-45 +20 



Table 1 

Upper bounds of the systematic error on the pr dependent single particle correction func- 
tion. Here Strack includes the uncertainties of the acceptance, dead areas, track matching 
cuts and the track reconstruction efficiency. The 5decay term accounts for the uncertainty of 
the decay correction. The effect of the momentum resolution contributes with dreso to the 
systematic error. Uncertainties due to potentially unsubtracted background are quantified 
by Sbgr- The total systematic error given in the last column is calculated as quadrature sum 
of the individual contributions. It is calculated separately for positive and negative errors. 



bin relative fraction 



N, 



part 



Ncoll N^ll'""'/^coll 2N^u/Npart 



1 

2 
3 
4 
5 



80-92% 
60-80% 
30-60% 
15-30% 
5-15% 
0-5% 



5.5± 2.6 
19.5± 3.5 
79± 4.6 
180± 6.6 
271±9 
348±10 



4.1±1.7 

20±6 
131±23 
406±46 
712±72 
1009±101 



246±98 
50.4±13 
7.68±1.1 
2.49±0.13 
1.41±0.03 
1 



1.5±0.5 
2.1 ±0.5 
3.4±0.6 
4.5±0.5 
5.2±0.6 
5.8±0.6 
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Table 2 

Number of participants and binary collisions and their systematic errors for the individual 

centrality selections used in this analysis. Also given is the ratio of the number of binary 
collisions for the most central sample relative to the one for each sample. The last column 
quantifies the ratio of binary colUsions to participant pairs. 
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